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ABSTRACT 

We apply our self-consistent accretion disk corona (ADC) model, with two different geometries, to the broad- 
band X-ray spectrum of the black hole candidate Cygnus X- 1 . As shown in a companion paper (Dove, Wilms, and 
Begelman), models where the Comptonizing medium is a slab surrounding the cold accretion disk cannot have a 
temperature higher than about 120 keV for optical depths greater than 0.2, resulting in spectra that are much softer 
than the observed 10-30keV spectrum of Cyg X-l. In addition, the slab geometry models predict a substantial 
"soft excess" at low energies, a feature not observed for Cyg X-l, and Fe Ka fluorescence lines that are stronger 
than observed. Previous Comptonization models in the literature invoke a slab geometry with the optical depth 
tt > 0.3 and the coronal temperature T c ~ 150 keV, but they are not self-consistent. Therefore, ADC models with 
a slab geometry are not appropriate for explaining the X-ray spectrum of Cyg X- 1 . Models with a spherical corona 
and an exterior disk, however, predict much higher self-consistent coronal temperatures than the slab geometry 
models. The higher coronal temperatures are due to the lower amount of reprocessing of coronal radiation in 
the accretion disk, giving rise to a lower Compton cooling rate. Therefore, for the sphere+disk geometry, the 
predicted spectrum can be hard enough to describe the observed X-ray continuum of Cyg X-l while predicting 
Fe fluorescence lines having an equivalent width of ~ 40 eV. Our best-fit parameter values for the sphere+disk 
geometry are tj sa 1.5 and T c « 90keV. 

Subject headings: radiation mechanisms: nonthermal - radiative transfer - X-rays: general - X-rays: binaries - 
accretion 



1. INTRODUCTION 

Due to its high X-ray brightness and its apparent large mass, 
Cygnus X-l is the most studied Galactic black-hole candidate 
(BHC). Although it is almost certain that the production of the 
high-energy radiation is associated with the accretion of matter, 
the exact details are uncertain. The source of accreting mat- 
ter is most probably a combination of Roche-lobe overflow and 
accretion out of the stellar wind of the compani on, the 09.7Iab 
supergiant HDE 226868 (pies & Bolton 1986|). Recent mass 



cal conditions in the Comptonizing medium), and he developed 
a more generalized theory. Using a composi te spectrum from 
EXOSAT, GRANAT, and OSSE observations, [Titarchuk (1994) 



determinations give a mass of about 1 8 M fl for the comp an- 
ion and 10 M for the comp a ct object (pterrero et al. 1995 , in 
agre ement with polan 199% |Sokolov 1987|; Bochkarev et al. 
1986; |Aab et al. 1984| ; |Hutchings 1978j and references therein), 
making the latter a good candidate for being a black hole. 

Cyg X-l is usually observed in the so-called "X-ray low, 7- 
ray high" state, which is roughly described by a power-law with 
a photon index a sa 1.7, modified by an e xponential cut-off 
with an e-fold energy of about lOOkeV (see |Oda 1977|; L iang 
& Nolan 1984; [Tanaka & Lewin 1995 for detailed reviews of 
the observations). The most popular model used for explaining 
the high-energy spectrum of Cyg X-l (and other BHCs), while 
in the "low state," involves inverse Comptonization of soft pho- 
tons by a hot plasma, usually r eferred to as a corona. While 
the Comptonization model by [Sunyaev & Titarchuk (1980) 



has been used widely i n the past (e.g., Grebenev et al. 1993 



Dobereiner et al. 1995[ ), [Titarchuk (1994)| showed that the typ 



ical fit-parameters found for Cyg X-l are outside the model's 
range of validity (and cannot be used for interpreting the physi- 



found that this new accretion disk corona (ADC) model, hav- 
ing an optical depth r e 1=3 0.6, a coronal temperature kT c 
150 keV, and a slab geometry allowed for a rough description 
of Cyg X-l. 

Although Comptonization apparently accounts for the gen- 
eral behavior of the observed X-ray spectrum of Cyg X-l, it 
appears that a reprocessing component is needed for the full 
description of the spectrum. Formal descriptions of the dif- 
ferent observations in terms of a reflection component with an 
underlying Comptonization spectrum give better fits to the data 
than do pure Comptonization mod els. The presence of a weak 
fluorescen ce line from neutral iron dKitamoto et al. 1990|; D one 
et al. 1992; Marshall et al. 1993| ; |Ebisawa et al. 1996| ) and the 
deviation of th e 2-60 keV- b and from a pure p ower law seen by 
HE AP 1-A2 (Unoue 1989|; |Done et al. 1992|; Gierliriski et al. 
1997), which is usually interpreted as a Compton reflection 
hump, both indicate the presence of cold or slightly ionized 
material present in the source, in addition to the Comptoniz- 
ing medium. Therefore, since the radiation processes of Comp- 
tonization and reprocessing of coronal radiation in the cold ac- 
cretion disk cover a large range of photon energies (< 2 keV to 
~ 250 keV) in the predicted X-ray spectrum, it is necessary to 
use-broad-band spectral observations to place meaningful con- 
straints on the models. 

As we discussed in our companion paper (Dove, Wilms, & 
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Begelman 1997, hereafter paper I), one serious drawback to 
solving the radiation transfer problem for an ADC with ana- 
lytic and "linear" Monte Carlo techniques (where each photon 
particle is propagated through a fixed background) is that the 
physical properties of the corona must be defined a priori. The 
coronal properties and the radiation field are strongly coupled 
through the processes of Compton cooling, Compton upscatter- 
ing, pair production and annihilation. Therefore, it is not clear 
whether the assumed coronal properties are self-consistent with 
the resulting radiation field. 

We have recently developed an ADC model in which the ra- 
diation field, the temperature and opacity of the corona, and 
the reprocessing of coronal radiation in the accretion disk are 
calculated self-consistently (paper I). We find that all previous 
ADC models having a slab geometry and used to describe the 
X-ray spectrum of Cyg X-l are not self-consistent, since the 
coronal temperatures and opacities are not within the physi- 
cally allowed ran ge of combinations of these parameters (see 
Fig. 2 of paper I). |Haardt et al. (1993~ claim that their "best-fit" 
model, having a coronal temperature kT c = 150keV and an 
optical depth Tt = 0.3, agree within the 3a confidence level 
of their allowed kT c — tj relation. However, due to the qual- 
ity of the data, the 3a contour covers a very large region of the 
kT c — tt parameter space, of which only a very small part is 
within the self-consistently allowed region. 

In this paper, we apply our self-consistent model to Cyg X- 
1, using a more compl ete composite spect rum than was avail- 
able to previous work (Wilms et al. 1997). We study two ge- 
ometries: (1) the standard slab geometry and (2) a spherical 
corona with an exterior accretion disk. At the time of writing, 
simultaneous X-ray spectra covering the range from 2 keV to 
250 keV were not generally available. We therefore had to use 
a non-simultaneous composite spectrum from several instru- 
ments to cover this energy range (each observation was taken 
while Cyg X-l was in its "low" state). A detailed description 
of the individual data comprising the composite spect rum and a 



preliminary analysis are given in Wilms et al. (1997) 



2. THE NUMERICAL CODE 

It appears impossible to produce an analytical model of an 
ADC that properly accounts for the self-consistent thermal and 
opacity distributions of the corona while solving correctly the 
radiative transfer problem of an angle-dependent, high-energy 
radiation field in a semi-relativistic, non-uniform plasma (in- 
cluding reprocessing of radiation in the accretion disk). There- 
fore, numerical methods are needed. For the computations pre- 
sented here, we use a non-linear Monte Carlo (NLMC) code 
to calculate self-consistently the temperature and opacity of the 
corona as well as the radiation field within and external to the 
corona. For the two geometries considered here, Compton scat- 
tering, photon-photon pair production, and pair annihilations 
are taken into account within the corona. For the reprocess- 
ing of radiation within the accretion disk, Compton scattering 
and photo-absorption, resulting in fluorescent line emission and 
thermal emission, are considered. We enforce charge neutrality 
by requiring n c — n + = n p , where n c , n + , and n p are the num- 
ber densities of electrons, positrons, and protons, respectively. 



For this paper, the spatial distribution of the seed opacity (the 
opacity not including the contribution due to electron-positron 
pairs) was assumed to be uniform. 

The free parameters of our models are: (1) the seed opacity 
of the corona, t p , (2) the temperature of the accretion disk as 
a function of radius, T BB (R), and (3) the compactness param- 
eter of the corona, l c . For a set of free parameters, the NLMC 
code is iterated until the temperature of the corona, T c , the total 
opacity, tj = r p + 2r + , and the radiation field have reached a 
steady state. Thermal equilibrium is achieved when the local 
heating rate is balanced by the Compton cooling rate. The pair 
opacity is determined by requiring pair-annihilations to balance 
pair-production. Once the system has reached steady state, the 
spectra of escaping radiation are recorded (in 10 bins of differ- 
ent inclination angles) until satisfactory statistics are obtained. 
A detailed description of the NLMC code and its use for ADC 
is given in paper I. 



3. SLAB GEOMETRY 

3.1. The Model 

The most common geometry used to model the spectra of 
BHCs is the slab (plane-parallel) geometry. For this geome- 
try, we assume that the accretion disk corona is situated above 
and below an optically thick, geometrically thin, cold accre- 
tion disk, that all corona and disk properties are constant with 
respect to radius, and that azimuthal symmetry applies. The 
coronal compactness parameter is defined to be 
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where zq is the scale-height of the corona, L c is the luminosity 
of the corona, i? c is the radius of the corona, i c dd is the Edding- 
ton luminosity, i? s is the Schwarzschild radius, ^ c is the rate of 
energy dissipation per unit area into the corona and h = zq /i? c . 
The specification of the physical mechanism heating the corona 
is not necessary. We simply assume that the energy is dissipated 
uniformly with respect to height. For a specific value of l c , the 
self-consist ent coronal prope rties are degenerate with respect to 
z and <f c ( [Stern et al. 1995[ ). 

As we discuss below, the predicted amount of thermal excess 
in the escaping spectrum of radiation strongly depends on the 
temperature of the accretion disk, T BB . A lower limit to the tem- 
perature can be estimated by assuming that the disk is heated 
solely by illumination from the corona. In this case, the tem- 
perature is approximated by balancing absorption with thermal 
emission: 



a sb T B 4 B = (1 - A)F C 



(2) 



where F c is the flux of coronal radiation incident onto the accre- 
tion disk (assumed here, by symmetry, to be equal to the flux of 
escaping radiation), cr s b is the Stefan-Boltzmann constant, and 
A is the angle-averaged albedo of the disk. Using equation ([[]), 

4 We note that, due to a bug in the treatment of reprocessing of coronal radiation, the conclusion made by us in Wilms et al. (1997) that the slab geometry model 
is able to describe the broad-band spectrum of Cyg X-l is erroneous. As we show below, the slab geometry model predicts a spectrum that is much softer than the 
observed spectrum. Due to the bug, the yields for Compton reflection/fluorescence and thermalization were reversed, causing < 20% (instead of > 80%) of the 
coronal radiation to be reprocessed into thermal radiation and leading to a porrp snonding errnnpons rnmptnn poo ling rate. As discussed in paper I, the current version 
of our code has been successfully compared to other reprocessing models (e.g., 3eorge, Nandra, & Fabian 199C ). 



3 



the disk temperature can be expressed as 
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We refer the reader to paper I for a more detailed discussion of 
the slab geometry models. 

3.2. Application to CygX-1 

As discussed in paper I, we find that the maximum self- 
consistent temperature of the corona is T MAX ~ 120 keV for 
total opacities tj > 0.2, regardless of the compactness pa- 
rameter l c and T BB . Once this maximum temperature, T MAX , 
is reached, higher heating rates result in a higher opacity due to 
pair production, causing more reprocessing of coronal radiation 
in the accretion disk. Since most of the reprocessed radiation 
is emitted as thermal radiation, the higher amount of reprocess- 
ing gives rise to higher Compton cooling rates. Therefore, once 
pair-production becomes important, an increase in the heating 
rate results in a decrease in the temperature. Fig. 2 of paper I 
shows the allowed parameter space of self-consistent corona 
temperatures and opacities. 

In Fig. |l] we show our by-eye "best-fit" of the predicted 
angle-averaged spectrum, modified by interstellar absorption, 
to the composite spectrum of Cyg X-l. For the slab geometry, 
the modeled spectrum shown is the hardest spectrum possible 
while having a cut-off energy E c > 100 keV. It is apparent, 
however, that the predicted spectrum is still much softer than 
the observed spectrum. In our grid of models, the range of seed 
opacities is 0.1 < r p < 2.0 and the range of coronal compact- 
ness parameters is 0.1 < l c < 10 3 (see paper I). Out of the en- 
tire grid, the hardest spectrum predicted by our self-consistent 
models has a photon index of a ~ 1.8, while Cyg X-l has a 
power -law index of ~ 1.6 — 1.7. As discussed in W ilms et al. 
(1997), The BBXRT archived data were deconvolved by assum- 
ing a power-law continuum, with a photon index of a = 1.62, 
absorbed by the interstellar medium with iVn = 6.0 x 10 21 
cm -2 . The TTM data were reduced by Bo rkous et al. (1995) 
the HEXE data by ^ 
by |Kurfess (1995)| 
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obereineret al. (1995)[ , and the OSSE data 
Since the models provided reasonable re- 
duced chi-squared values (< 1.5), the unfolded data should be 
at least a good approximation to the actual spectrum. There- 
fore, even though we are comparing our predicted spectra with 
unfolded data, which are dependent on the assumed model, the 
result that the predicted spectrum is always softer than the ob- 
served spectrum is not sensitive to the deconvolution procedure. 

According to our linear (and therefore not self-consistent) 
MC simulations, simulated spectra with a < 1.7, while T c > 
100 keV, are achieved only for models with r > 0.3. These 
values are consistent with past work, in which the spectrum of 
Cyg X-l has been described by ADC mod els with a slab geom- 



etry (Haardt et al. 1993; ritarchuk 1994). Such models, how 



ever, are not self-consistent, for the temperature and opacity 
values lie within the forbidden region of Fig. 2 of paper I. We 
remind the reader that the predicted maximum temperatures are 
indeed upper limits, for our models do not include additional 



cooling mechanisms such as bremsstrahlung. (This mechanism 
was found to be negligible for the models of interest, but its 
inclusion would reduce the maximum temperatures allowed.) 

In addition to the simulated spectra having too soft a power- 
law component, our models predict a very large Fe Ka equiva- 
lent width (EW). Our "best-fit" slab model predicts an EW ~ 
120 eV , considerably larger than the observed value of EW < 
30 eV ( Ebisawa et al. 199q ). Finally, the slab models also pre- 
dict an excess of radiation for energies E < IkeV (i.e., a 
"soft excess"). Since the "best-fit" slab models are optically 
thin, most of the thermal emission emitted by the disk escapes 
the system without interacting with the corona. We have ac- 
counted for Galactic absorption by using th e observed colum n 
density of hydrogen, iVn = 6 x 10 21 cm~ 2 (Wilms et al. 1997). 



Clearly, this is not sufficient to "hide" the soft excess for the 
200 eV disk models. Models with a higher value of iVn predict 
too low of a flux for energies below 1 keV, so this soft-excess 
problem cannot be solved by having having arbitrarily large hy- 
drogen column densities in the model. We note, however, that 
the BBXR T observation is uniq ue since a soft-excess was not 
observed (Marshall et al. 1993 ). Since Cyg X-l usually has 



a soft-excess component (e.g., Balucihska & Hasinger 1991 



Balucinska-Church et al. 1995; Ebisawa et al. 1996), the dis 



crepancy between the predicted flux and the observed flux at 
low energies should not be considered a serious problem by it- 
self. However, it does appear that our "best-fit" slab geometry 
model predicts a higher soft-excess component than observed 
by the AS CA GIS detector for November 11, 1993 (E bisawa 
et al. 1996). 

We conclude that self- consistent accretion disk corona mod- 
els with a slab geometry are not capable of reproducing the 
observed broad X-ray spectrum of Cyg X-l. 

4. SPHERE+DISK GEOMETRY 

4.1. The Model 

The main limitation of the slab geometry models, in regard 
to explaining the high-energy spectra of Cyg X-l, is that there 
is too much reprocessing of coronal radiation in the cold ac- 
cretion disk. Therefore, models with the reprocessing matter 
having a smaller covering fraction (as seen from the corona) 
are more likely to explain the observations. One such geom- 
etry is a combination of a spherical corona with an optically 
thick, geometrically thin, cold accretion disk. Here, the accre- 
tion disk is assumed to be exterior to the corona, as shown in 
Fig. ^. This geometry is very similar to the two-te mperature 
accretion disk model of Shapiro & Lightman (1976), although, 
in our model, the proton temperature is assumed to be equal 
to the electron temperature. This geometry is also similar to 
the adve ction-dominated disk models by Narayan an d collab- 
orato rs ([Narayan & Yi 19954 Narayan & Yi 1995b|; N arayan 
1996, and references therein). However, in our models, the seed 
photons for Comptonization are produced by the accretion disk 
through thermal emission rather than bremsstrahlung and syn- 
chotron radiation within the corona. The accretion disk extends 
from R c to i?d where R c is the radius of the sphere and i?d is 
the outer radius of the disk. We define the ratio of the outer disk 
radius to the corona radius as a = i?d/i? c . For this paper, we 
set a = 10, but the main results are not sensitive to this value 



(see § X.2). 



4.1.1. Energetics 
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Fig. 1 . — Slab geometry: Comparison of the predicted spectrum of escaping radiation, modified by interstellar absorption, to the 
composite X-ray spectrum of Cyg X- 1 . For our model, T BB = 200eV, r T = 0.28, T c = llOkeV, and N H = 6 x 10 21 cm~ 2 . 




Fig. 2. — The sphere+disk geometry. Spherical corona with exterior accretion disk. The corona has a radius R c , and the accretion 
disk (assumed to be optically thick) has an inner radius R c and an outer radius i?d- The corona has been divided into several zones, 
such that non-uniform temperature and opacity structures can be studied. For this paper, we assume uniform heating with volume 
within the corona. 
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We define the coronal compactness parameter by 



m e c 3 Rc 
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Similar to the slab geometry models, we allow for intrinsic 
emission of thermal radiation by the accretion disk (in addi- 
tion to emission due to the reprocessing of coronal radiation). 
We define fP G to be the rate at which gravitational energy is 
dissipated directly into the corona, where Pq is the total rate of 
gravitational energy dissipated into the system. Consequently, 
(1 — /)Pg is the rate °f energy dissipation into the disk. The 
total luminosity of the accretion disk is given by 



L d = (1 - f)P G + L 



absi 



(5) 



where L a b s is the rate at which energy is absorbed (not re- 
flected) by the disk due to the reprocessing of radiation emit- 
ted by the corona. Equation (|j) can be expressed in terms of 
compactness parameters as 



k = (1 - f)l G + I 



abs ; 



(6) 



where l d = (a T /m e c 3 )L d /R c , l G = (a T /m e c 3 )PG/R c , and 
'abs = (o"t /m e c 3 ) L a bs / Rc- As with the slab geometry models, 
we constrain the models by setting (1 — f)l G = 1 for all of the 
simulations. With this choice, / is given by 



lc 



1 + lc 



(7) 



As discussed in paper I, setting (1 — f)l G to unity allows us to 
consider models where 0.01 < / < 1.0, but models with other 
values of (1 — f)l G yield the same ranges of self-consistent 
coronal temperatures and opacities. 

4.1.2. Disk Temperature 

As in the slab-geometry model, we can estimate a lower limit 
to the disk temperature by assuming that the disk is heated 
solely by illumination from the corona; 



7ri? c 2 (a 2 - l)a sb (T BB ) 4 = / d ( )(l - A) 



(8) 



where L c is the total luminosity leaving the corona and /d(») is 
the fraction of photons leaving the corona that is reprocessed in 
the disk. This fraction can be approximated by assuming that 
the photons leaving the corona are uniformly distributed over 
the surface of the corona with an isotropic distribution. With 
these assumptions, the fraction of escaping photons that hits 
the disk is given by 
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where (Aft) is the average solid angle of the disk as seen from 
the surface of the corona, 
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m = arcsin(l/a), and is the angle between the normal of 
the coronal surface and the plane of the disk. In Fig. ||, we 
show how the covering fraction depends on a. It is apparent that 
fd(a) is nearly constant for a > 10, and therefore the amount 
of reprocessing of radiation within the accretion disk is insen- 
sitive to a. The maximum value of f d is « 1/ 3. In contrast, 
for a slab geometry, all of the downward directed radiation at 
the base of the corona interacts with the accretion disk. There- 
fore, the slab geometry models predict much more prominent 
reprocessing features in the spectrum of escaping radiation. 

For large accretion disks (i.e., a 3> 1), the local disk tem- 
perature can vary by more than an order of magnitude between 
the inner and outer radius. Rather than using an average disk 
temperature for calculating the spectrum emitted by the disk, a 
more proper treatment is to take into account the temperature 
structure. We can estimate the disk temperature as a function of 
radius by equating the flux of absorbed coronal radiation with 
thermal emission, 



a sb T B 4 B (r) 2nrdr = -(1 - A)L c df d , 



(10) 



where r = R/R c , df d (r)/dr is the differential covering frac- 
tion of a ring of radius r, and f d (r) is determined numerically 
using equation (^). Using equation (Q), the temperature can fi- 
nally be expressed as 



T BB (r) = 150 
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In Fig. [i]we show how T BB (r) varies with r. From a numerical 
fit, we find that T BB (r) roughly decreases as T BB (r) oc r 
as compared to a po wer-law of 3/4 correspondi ng to the stan- 
dard a-disk models ( [Shakura & Sunyaev 1973[ ). It should be 
remembered, however, that if any accretion-energy were dis- 
sipated directly into the accretion disk (/ < 1), then the disk 
temperature would be higher than the values given above. To 
address this possibility, we consider a simple model where, for 
R > R c , the accretion d isk behaves as the standard a disk 
(Shakura & Sunyaev 1973), while, for R < R c , all of the accre- 



tion energy is dissipated into the corona. In order to determine 
the relative importance of the internal dissipation of energy, 
compared to the contribution from coronal illumination, we es- 
timate the disk temperature by neglecting the contribution from 
illumination. In addition, we make the simplifying assumption 
that the disk is optically thick and radiates as a blackbody such 
that the emitted flux is equal to the dissipation rate, 



a sb T£ B (R) = D(R) 



( 3GMM 
1 8irR 3 



1/4 



(12) 



where we have assumed R^> R s (Frank, King, & Raine, 1992). 
Defining L acc = r]Mc 2 , where r\ is the efficiency by which the 
accretion process converts gravitational energy into radiation, 
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Fig. 3. — The covering fraction, /d(a), of the accretion disk with an outer radius a (normalized to the coronal radius), averaged over 
the surface of the corona (spherical geometry). 
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Fig. 4. — The temperature of the accretion disk as a function of radius. The temperature is calculated by assuming the disk is heated 
by illumination of coronal radiation and cooled by thermal emission (equation (Mb). Here, l c = 25, M = 1OM , and R c = 100i? s . 
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and using the standard definition of the Eddington luminosity, 
Zvcdd> we can express the disk temperature as 



kT BB (R) = k 
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By comparing equation ( |13| ) with equation ( |1 1| ) (or Fig. |J), it 
appears that the heating rate due to the direct dissipation of 
accretion energy is comparable to that from coronal illumi- 
nation. A proper treatment of the disk temperature structure 
would be to equate thermal emission with both heating mecha- 
nisms. Such a treatment is outside the scope of this paper; with 
a heating rate equal to twice the value used in equation (|T3|), 
the disk temperature will be roughly 1.2 times higher than the 
value given in equation ( |l3| ) (T BB is proportional to the total 
heating rate per unit area). Therefore, the inclusion of direct 
energy dissipation into the disk changes the disk's total lumi- 
nosity, but does not significantly change the temperature profile 
of the disk. 

For simplicity, we have assumed that the intrinsic flux of 
thermal radiation is constant with respect to radius (the emis- 
sion due to reprocessing, however, is determined locally by 
equating emission with absorption). This approximation is not 
too important since the emission due to reprocessing is higher 
than the emission due to internal dissipation for the models 
considered in this paper. In regards to the soft excess issue, 
for a ss 10, these models slightly overestimate the predicted 
amount of soft excess by less than 10%. 

4. 1 .3. Covering Fraction of Corona 

In all of our models, the seed-photons are produced within 
the cold accretion disk (we neglect bremsstrahlung radiation). 
For this geometry, however, not all seed photons have to prop- 
agate through the corona prior to escaping the system. There- 
fore, even for optically thick models, a very large soft excess 
of radiation is always predicted. The covering fraction of the 
corona, evaluated on the accretion disk at a distance d from the 
center, is given by 
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where R c is the radius of the corona. As shown in Fig. the 
covering fraction of the corona rapidly decreases with increas- 
ing R. In order to estimate the fraction of the thermal radiation 
emitted by the accretion disk that interacts with the corona, we 
average f c over the accretion disk, weighted by the thermal flux, 
F(r), 



(M(a) 



f" 2nrf c (r)F(r)dr 
2nrF(r)dr 



(15) 



For an accretion disk heated solely by illumination, we found 
that T(r) oc r -11 (Fig. Q). On the other hand, for the case 
where the disk is heated solely by internal dissipation, i.e., the 
standard a-disk, T{r) oc r -3 / 4 (Shakura & Sunyaev 1973). 
For a rough estimate of (/ c ) when both disk heating mecha- 
nisms are important, we assume F(r) = kT 4 (r) oc r~ 35 . In 



Fig. ^| we show how this averaged covering fraction depends on 
a. For a > 2, (/ c ) < 0.2, and at least 80% of the thermal ra- 
diation escapes the system without interacting with the corona, 
even if the corona is optically thick. In contrast, for the slab 
geometry models, all of the thermal radiation must propagate 
through the corona prior to escaping the system. Therefore, 
the sphere+disk models always predict that a large amount of 
thermal radiation, relative to the amount of Comptonized radi- 
ation, will escape the system. Unless this thermal radiation is 
absorbed by the interstellar matter (ISM), it will appear as a soft 
excess of radiation. 

For Tt > 1, the self-consistent temperature structure varies 
by roughly 30% for the case where the heating rate is uniform 
with volume. For these optical depths, the outer shells of the 
corona are the hottest while the center region is the coldest. 
This variation is not large enough such that the predicted spec- 
tra contain a "hardenin g" feature that mimics t he hardening due 
to Compton reflection (Skibo & Dermer 1995). 

For the sphere+disk geometry, the accretion disk receives 
less illumination from the corona since the covering fraction is 
lower than in the slab geometry case, while the surface area of 



the disk is larger. In § X. 1 , we estimated the temperature struc- 
ture due to either coronal illumination or internal dissipation of 
accretion energy (equation n~3t). For either case, 



T BB (r) oc r~ 



(16) 



where r = R/R s , J f» 1.1 for heating by coronal illumination 
and 7 = 0.75 for the standard a-disk. The flux of radiation 
emitted by an accretion disk is approximated by 



F(E) 



B E (T)rdr, 



(17) 



where B E (T) is the Planckian distribution. In Fig. ^, we show 
the corresponding spectrum, F(E), for 7 = 1.1, and we com- 
pare F(E) to a Planckian distribution corresponding to the tem- 
peratures (T BB ) = 60 eV and 80 eV. It is seen that the Planckian 
distributions approximate the high-energy tail fairly- well. The 
disagreement at energies lower than 100 eV is irrelevant since 
this radiation is efficiently absorbed by the ISM. In addition, 
since the temperature decreases rapidly with radius, most of the 
radiation having an energy E > 100 eV is emitted within the in- 
ner region of the disk, and the exact value of 7 is unimportant. 
Therefore, since the Comptonized spectrum is insensitive to the 
exact shape of the spectrum of seed photons, we used a Planck- 
ian distribution corresponding to a single disk temperature of 
80 eV to approximate the exact integrated disk spectrum. 
4.2. Application to Cyg X-l 

We computed a grid of models, in which the range of seed 
opacities is 0.5 < t p < 4.0 and the range of coronal compact- 
ness parameters is 0.1 < l c < 100. In Fig. ^, for kT BB — 50eV, 
we compare the predicted spectrum of escaping radiation, mod- 
ified by interstellar absorption, of our "best-fit" model to the 
composite spectrum of Cyg X-l. Since the seed photons enter 
the corona from the exterior of the sphere, the effective opti- 
cal depth of the sphere is reduced. Therefore, the total optical 
depth (as defined by tj = R c crT;(n e + n p ), where n e = n p + n + ) 
must be higher than that for the slab geometries in order for the 
models to predict similar power-laws. Our "best-fit" model has 
a total optical depth tj = 1.5 and an average coronal tempera- 
ture (T c ) = 90keV. 
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Fig. 5. — The covering fraction of the spherical corona (with radius i? c ) evaluated on the disk at a radius R. 
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Fig. 6. — The covering fraction of the spherical corona, averaged over the accretion disk with outer radius a = R&/R c and weighted 
by the flux of thermal radiation, F(r) oc r~ 3 - 5 . 
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Fig. 7. — The spectrum of thermal radiation emitted from an accretion disk with a radial temperature structure, T BB (r), as given by 
equation [ll] (solid line). Here, l c = 23, R c = 100i? s , and M = 7M Q . Dashed line is a Planckian distribution corresponding to 
a temperature T BB — 60 eV. Dashed-dot-dashed line is a Planckian distribution corresponding to a temperature T BB — 80 eV. All 
distributions have been normalized to unity. 
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Fig. 8. — Comparison of the predicted spectrum of escaping radiation, modified by interstellar absorption, with the X-ray spectrum 
of Cyg X- 1 . This is for the sphere+disk geometry. Here, T BB = 50 eV, t t = 1 .5, T c = 90 keV, and iV H = 6 x 10 21 cm" 2 . 
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Fig. 9. — Allowed temperature and opacity regime for self-consistent ADC models with a sphere+disk geometry. Solid line is 
derived from a "fit by eye" to the numerical results. For a given total optical depth, temperatures above the solid line are not possible. 
For all models, the blackbody temperature of the disk is kT BB = 50 eV. Different symbols represent models with different seed 
opacities, while different points having the same symbol represent different coronal compactness parameters. 



With this geometry, the corona is able to reach much higher 
temperatures as compared to the slab models because the 
corona is "photon-starved," i.e., the luminosity of the seed 
photons is much less than the luminosity of the hard X-rays 
(Zdziarski et al. 1990). As discussed above, most photons that 
are reprocessed within the accretion disk do not re-enter the 
corona. Therefore, the Compton cooling mechanism that pro- 
hibits slab geometry coronae from having high temperatures 
is not as efficient in keeping the spherical coronae cool. In 
Fig. !| we show the allowed regime of optical depths and aver- 
age corona temperatures for the sphere+disk geometry. Since 
the coronae are allowed to be much hotter than the slab ge- 
ometry models, there are self-consistent {T c ,tj} combinations 
such that the observed power-law and cut-off in BHCs can be 
reproduced. Since more than 90% of the thermal radiation es- 
capes the system without passing through the corona, the soft 
excess is very large for these models. As in the case with the 
slab-geometry models, Galactic absorption is not sufficient to 
"hide" the soft excess for the 200 eV disk models. In Fig. ^[ 
for kT BB = 50 eV, we compare our best-fit model to the spec- 
trum of Cyg X-l. Here, the predicted thermal excess appears to 
be consistent with the observational data, and the entire-broad 
band spectrum is adequately described by the model. As dis- 
cussed in the Appendix, even though we are fitting 'unfolded' 
data, we are confident that this model provides a good represen- 
tation of the intrinsic photon-spectrum of Cyg X-l. However, 
the main point of this paper is to motivate a more detailed study 
of the sphere+disk model, using a more rigorous method of data 
analysis. 

A nice feature of the sphere+disk model is that it naturally 
predicts an "effective" disk temperature to be much lower than 
in the slab geometry models. In fact, due to the lower disk 
temperature, this model appears to predict a soft excess that is 



consistent with the observations (Baluciriska & Hasinger 1991 
Balucihska-Church et al. 1995). The quality of the observa 



tions, however, is not too good for E < 2keV, which is un- 
fortunate since this is the energy range where the spectru m is 
most sensitive to the disk temperature. Also, as noted in § 3.2, 
Cyg X-l usually contains a soft-excess component. We plan 
on comparing this model to better observational data in this en- 
ergy range (e.g., RXTE and ASCA data) in order to constrain 
this model further. 

To test whether the predicted thermal radiation from the ac- 
cretion disk is observable in the UV, we have compared the UV 
flux predicted by our model for the spectral band from 1200 A 
to 1950 A with the flux observed in this band by IUE. Our pre- 
dicted flux is two to three orders of magnitude smaller than the 
flux in an IUE Low Dispersion spectrum of HDE 226868 made 
in June 1980. This is consistent with earlier findings that the 
UV spectrum of the system HDE 2 26868 — Cyg X-l is dom- 
inated by e mission from the O star ([Treves et al. 1980|; Pravdo 
et al. 1980). Therefore, IUE observations cannot be used to 
constrain the disk temperature of our models. 

Due to the smaller covering fraction of the accretion disk, 
the models predict weaker reprocessing features in the escap- 
ing radiation field as compared to the slab models. In principle, 
in modeling BHCs, the value of a can be constrained by the 
equivalent width of the Fe Ka fluorescent line, which is pro- 
portional to the solid angle of the disk, /d (a) . For Cyg X- 1 , the 
EW of the Fe Ka line is measured to be EW < 70 eV (E bisawa 
et al. 1996; pierlihski et al. 1997[ ). With a =10 (correspond- 
ing to a disk covering fraction /<j w 0.3), our best fit models 
predict an EW of ~ 60 eV. Therefore, models with smaller /<j 
will predict EWs that are too low to be consistent with the ob- 
servations. For /d > 0.25, a > 3 (Fig. |J). On the other hand, 
since /d(a) does not vary significantly for a > 10, this method 
cannot provide an upper limit on the size of the accretion disk. 

5. DISCUSSION 

Our self-consistent slab accretion disk corona models are un- 
able to explain the broad-band X-ray spectra of Cyg X-l. The 
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modeled coronae are either too cold or have an optical depth 
that is too small, resulting in spectra that are much softer than 
the observed spectrum. All previous ADC models that have 
successfully described the spectra of Cyg X-l have tempera- 
ture and opacity values that are outside the allowed region of 
self-consistent values. In addition, as discussed in paper I, the 
predicted angle-averaged equivalent width (EW) of the Fe K a 
fluorescent line is EW > 150 eV, a value roughly three times 
high er than the measured value for Cyg X-l (E bisawa et al. 
1996), indicating that the covering fraction of the reprocessing 
material is il/2n ~ 0.3. 

We believe that there is no way around these shortcomings of 
the slab geometry model and that these ADC models are not the 
appropriate models for explaining the high-energy spectra of 
BHCs. This claim is in agreement with Gierliriski et al. (1997), 
who, by fitting the joint Ginga-OSSE observation of Cyg X-l 
with a power law plus a reflection component, find that the solid 
angle of the cold medium is ~ 0.3 x 2tt, significantly lower than 
the a ngle corresponding to t h e slab geometry, a result consistent 
with Ebisawa et al. (1996) Gierlinski et al. (1997) also argue 
that the corona must be "soft-photon starved," which rules out 
a slab geometry. 

For ADC models with a sphere+disk geometry, there are 
self-consistent temperature and opacity combinations such that 
the power-law and cut-off portions of the spectra of BHCs can 
be explained. In addition, due to the small covering fraction of 
the disk as seen from the corona, it is possible that the "effec- 
tive" disk temperature (a flux-weighted average disk temper- 
ature) can be as low as < 80 eV. With such low values, the 
models predict a thermal excess (or lack thereof) that appears 
to be consistent with the observations of Cyg X-l since Galac- 
tic absorption is more efficient at these lower energies. It is still 
unclear, however, whether an ADC model havin g this geome- 
try is stable, as the Shapiro & Lightman (1976) model suffers 
from a thermal instability ( firan 1978 ). In our models, we can 
only assume a corona/disk morphology, but we cannot deter- 
mine whether this morphology is thermally or dynamically sta- 
ble. The sphere+disk geometry does approximate the geometry 



cently implemented the sphere+disk geometry spectral grid a s 
an user-defined model for XSPEC (version 9, Arnaud 1996), 



corres ponding to the advection-dominated models of Marayan 



(1996, and references therein) which appear to be thermally sta- 
ble. In addition, our best-fit models for Cyg X-l have an optical 
depth Tt ~ 1.5 and a coronal temperature kT c w 90keV, val- 
ues that are consistent with the advection-dominated models if 
the mass accretion rate M ~ 0.1 — 0.03Mc dd. where M e dd is 
the Eddington accretion rate ( Narayan 1996| ). 

The results obtained here for Cyg X- 1 are also applicable to 
other low state BHCs. In Fig. [l(] we compare the spectra of the 
black hole X-ray transients GS2023+338 and GRO J0422+32 
to the predicted spectra from our sphere+disk geometry model. 
Occasionally black hole X-ray transients have spectra that are 
dominated by a strong soft-excess at lower energies and a pure 
power-law component at high energies. In the low state, how- 
ever, GS2023+338 (=V404 Cyg) and GRO J0422+32 did not 
contain a soft excess, and their spectra of are very similar to 
that of the low-state of Cyg X-l. Therefore, it appears that 
photon starved ADC models, such as the sphere+disk geome- 
try, are applicable in explaining the observed X-ray spectra of 
many BHCs. 

Our "fits" to the unfolded data should be interpreted as a mo- 
tivation to consider this model in much more detail. For exam- 
ple, we have recently implemented our sphere+disk model in 
tabular form for the XSPEC user-defined model. We have re- 



and are currently analyzing new data of Cyg X-l and other 
BHCs from the RXTE satellite. Our preliminary results are in 
agreement with the results presented in this paper. 
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APPENDIX A 

Comparisons of our models with observations have relied on 
the photon-spectrum presented in Fig. ^. Although the use of 
"unfolded" X-ray spectra is very common in the astronomi- 
cal literature, it is important to keep in mind how these spec- 
tra are prepared. Due to the poor energy resolution of today's 
X-ray detectors, the inv ersion of the detector r esponse-matrix 
is generally impossible (Blisset & Cruise 1979, and references 
therein). Thus, in order to get an estimate of the intrinsic pho- 
ton spectrum from a source, a spectral model has to be fit to 
the spectrum observed by the X-ray detector, usually using a 
X 2 -reduction method and an estimate for the detector response 
matrix (Gorenstein, Gursky, & Garmire, 1968). If the x 2 -value 
resulting from the fitting process is small, then it is assumed 
that the modeled photon-spectrum is a good representation of 
the observed photon-spectrum. In order to produce an "un- 
folded" spectrum, the additional assumption is made that the 
residual (for a given detector energy channel) between the mod- 
eled spectrum (folded through the response matrix) and the ob- 
served count-rate is a good representation of the deviation be- 
tween the spectral model and the intrinsic photon-spectrum (in 
energy space). Thus, at a given energy, the product of the mod- 
eled photon-flux and the ratio between the observed count-rate 
and the (estimated) model count-rate is what is called the "un- 
folded" spectrum. 

Although a small x 2 -value of the fit to the "unfolded" spec- 
trum usually implies that the model is an adequate description 
of the intrinsic spectrum, sometimes this is not true. Fitting our 
models to the observed data in detector-space (using the method 
described above) would, in principle, give the best method of 
determining the "goodness of fit." However, the CPU time 
needed to obtain model spectra with a signal to noise ratio that 
is suitable for spectral fitting is prohibitively high, and we were 
forced to use a simplified approach; for a particular detector, 
we simulated an observation using our modeled spectrum as the 
intrinsic spectrum, the correct detector response matrix, back- 
ground estimate, and observational Poisson noise. This sim- 
ulated spectrum was then fit with several spectral models that 
are typically used to describe the X-ray spectra of BHCs (e.g., 
a power-law with an exponential cutoff). 

We performed these tests using several typical X-ray detec- 
tors (ROSAT, XTE-PCA, HEXE) and using typical count-rates 
and observation times for these detectors. In all of the tests the 
fit-parameters found using our model spectrum were in agree- 
ment with typical fit-parameters found in the real observation 
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Fig. 10. — Comparison s between the composite spectra of GRO J0422+32 ( Kroeger et al. 1996) and V404 Cyg ( in' t Zand et al. 1993 



Dobereiner et al. 1994; June 20, 1989 spectrum) and the predicted spectra from the best-fit sphere+disk models. The spectra have 
been scaled arbitrarily. For V404 Cyg, the parameter values of our best-fit model are kT BB = 52 keV, tt = 3.4, and iVn = 2.4 x 10 22 
cm 2 . For GRO J0422+32, kT BB = 75keV, r T = 1.9, and N H = 1.6 x 10 22 cm 2 . 



of Cyg X-l. In the simulated XTE-PCA observation, e.g., the 
model spectrum below 30keV can be roughly described by an 
absorbed power-law with an index o f a pa 1.62. Th is in- 



dex is similar to those f ound, e.g., by Done et al. (1992) and 



Dobereiner et al. (1995) In the residuals of the simulation the 



Fe line and the reflection hump are clearly visible. Using a 
more appropriate model that incorporates reflection, the simu- 
lated spectrum is best described by reflection of a power-law of 
index T w 1.72 off a cold accretion disk with covering factor 



0.4. Using an assumed length of 5 ksec for the ROSAT-PSPC 
observation it is possible to see evidence for the soft-excess in 
the simulated observation. According to this simulated obser- 
vation, the soft-excess as seen by ROSAT PSPC can be de- 
scribed by a black-body with a temperature of about 100 eV 
or less, which is roughly compar able with the observations by 
Baluciriska-Church et al. (1995) . Therefore, we are confident 
that the model-spectrum of Fig. || is a good representation of 
the intrinsic photon-spectrum of Cyg X-l. 
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